DETERMINATION OF \V US \ FROM HADRONIC r DECAYS 
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The recent update of the strange spectral function and the moments of the invariant mass 
distribution by the OPAL collaboration from hadronic r decay data are employed to determine 
\V u s\ as well as m s . Our result, \V US \ = 0.2208 ± 0.0034, is competitive to the standard 
extraction of \V U3 \ from K e s decays and to the new proposals to determine it. Furthermore, 
the error associated to our determination of \V U3 \ can be reduced in the future since it is 
dominated by the experimental uncertainty that will be eventually much improved by the B- 
factories hadronic r data. Another improvement that can be performed is the simultaneous fit 
of both \V US I and m 3 to a set of moments of the hadronic r decays invariant mass distribution, 
which will provide even a more accurate determination of both parameters. 



1 Introduction 

Already in the past, hadronic r decays have served as an interesting source to study low energy 
QCD under rather clean condi tion s^ and obtain information on parameters of the Standard 
Model, like the strong coupling^^, the strange quark mass or some non-perturbative conden- 
sates. The determination of a s , for example, has been performed using r decay data with a 
precision competitive to the current world average. 
At present, the hadronic r decay width 

_ T [r~ -> hadrons(7)] 

T ~ r[7--e-F e i>T(7)] 1 ' 

as well as invariant mass distributions, have reached a high precision status thanks to the data 
of the LEP experiments ALEPH^ and OPAL^ at CERN and the CESR experiment CLEO El 
at Cornell. 

In particular, the experimental measurements of the strange spectral function l ^ l ^ l ^ l have 
allowed the analysis of the SU(3) breaking corrections in the semi-inclusive r decay width into 



Cabibbo suppressed modes with strange particles, providing a way to compute the strange 
quark mass m s (M T ) l o l ^ l 1 ^ 11 ! 1 ^ 1 ^ ! 14 ! 10 !. Recently, we have pointed out that this approach to 
obtain the strange mass from the hadronic r decay s de pends sensitively on the modulus of the 
Cabibbo-Kobayashi-Maskawa matrix element |V^ s |^l It appears then natural to turn things 
around and, with an input for m s obtained from other sources, to actually determine \ V US \. The 
great advantage of the determination of this CKM matrix element from r decays in comparison 
with other calculations is that the experimental uncertainty, which as we will see is the main 
source of error in the calculation, is expected to be reduced drastically at the present B-factories 
BABAR and BELLE. 



2 Theoretical framework 

The main quantity of interest for the following analysis is the hadronic decay rate of the r 
lepton defined in The basic objects one needs to perform the QCD analysis of R T and 
related observables are Green's two-point functions for vector V/j = qff^qj and axial-vector 
A ij = 9i7 M 759j color singlets, 



(2) 



H%(?) =if d 4 xe^(0|r([^]t(x)^(0)) |0), 

:|d 4 xe^(o|r(K/(x)^(o)) |o). 



The subscripts i, j denote light quark flavors (up, down and strange). These correlators admit 
the Lorentz decompositions 



(q) 



9^q Z + q»q v ) ^, V/A tf) + ?V %y t A^ 



(3) 



where the superscripts in the transverse and longitudinal components denote the spin J = 1 (T) 
and J = (L) in the hadronic rest frame. Theoretically, R T can be expressed as an integral of 
the imaginary part of these correlators over the invariant mass s = p 2 of the final state hadrons^ 







2 r / 2s \ 

1 + ^2 )Imn T (s)+Imn L (s) 



where the appropriate combinations of two-point correlation functions are 



tA,J, 



+ Wu. 



(4) 



(5) 



with Vij being the corresponding matrix elements of the CKM matrix. 

Experimentally, one can disentangle vector from axialvector contributions in the Cabibbo- 
allowed (ud) sector, whereas such a separation is problematic in the Cabibbo-suppressed (us) 
sector. We can then decompose R T both experimentally and theoretically into 



R T = R T ,v + Rt,a + Rt,s- 



(6) 



Additional information can be inferred from the measured invariant mass distribution of the 
final state hadrons, through the analysis of the moments 



R (k,l) 



Mi 



ds 1 



\ k f s V dR T 



M?J \M^J ds 



(7) 



that can be calculated in analogy to the r decay rate (R T = R T ) and, in particular, can be 
also decomposed into Cabibbo-allowed and Cabibbo suppressed contributions. 



The theoretical study of R T and its moments is based on the Operator Product Expansion 
(OPE) of the relevant correlators. In this framework, the moments R^' 1 ^ can be written as 



4*'') = N c S EW {(\V ud \ 2 + \V u , 12 > 



) [l + 5^} + £ [\V ud \ 2 & iD) + \V US \ 2 &^] } . (8) 



D>2 



The electroweak radiative correction Sew = 1-0201 ± 0.0003 ^ has been pulled out explicitly 
and <5( fc '^(°) denotes the purely perturbative dimension-zero contribution, that is the same for the 

(k l)(D) 

Cabibbo-allowed and Cabibbo-suppressed contributions. The symbols o^-' stand for higher 
dimensional corrections i n the O PE from dimension D > 2 operators which contain implicit 
l/M^P suppression factors P M 1 1 1 1 7 1 The most important corrections are the dimension D = 2 
proportional to m 2 . and the dimension D = 4 proportional to m s (qq). 

The sep arate measurement of strange and non-strange contributions to the decay width of 
the r leptod^EI allows one to pin down the flavour SU(3)-breaking effects, dominantly induced 
by the strange quark mass, through the differences 



SR^' l) 



R 



(A,0 
t,V+A 



R 



(k,i) 

T,S 



W I 2 

I ' us I 



N, 



Sew £ [& (D) ~ 6% im 



D>2 



0) 



Many theoretical uncertainties drop out in these observables since they vanish in the SU(3) 
limit. 



3 Calculation of \V US \ 

The large sensitivity of the SU(3) breaking quantities 6R+ allows us to obtain a determination 
of the CKM matrix element \V US \ using as input a fixed value of m s . Since the sensitivity to 
I'Ktal is strongest for the moment with k = I = 0, where also the theoretical uncertainties are 
smallest, we used this moment in our calculation. From Q, 



R (o,o) 

l ry+A _ $ R (0,0) 



In the OPE of MjJti? 

in © we include the dimension two corrections S^' 1 ^^ that are known 
at 0(a») for J = L component and at Oia 2 ) for J = L + T component -see™ for references. 
The O{o?) J = L + T are also known^^but the results obtained in^3 and discussed here do 
not contain the se corrections. We leave a full analysis which take them into account for a future 
publication ^9. Dimension four corrections are fully included while the dimension six 

corrections were estimated to be of the order or smaller than the error of the dimension 

four El. 

An extensive analysis of the perturbative series for the dimension two corrections was done 
in 1*3 The conclusions there were that while the perturbative series for J = L + T converges very 
well the one for the J = L behaves very badly, adding important uncertainties to the theoretical 

(k I) 

calculation of 5R)- .A natural remedy to solve this problem is to replace the QCD expressions 
of scalar and pseudoscalar correlators by corresponding phenomenological hadronic parametriza- 
tions, much more precise than their QCD counterpart'^ due to the fact that it is dominated by 
far by the well known kaon pole -see^^ for details. The results we obtained within the QCD 
correlators and within the phenomenological hadronic parametrized correlators are showed in 
Table 1. The longitudinal contributions calculated with the two different descriptions of the 



Table 1: Comparison between the OPE and the phenomenological hadronic parametrizations explained in the 

text for the longitudinal component of R} !?,,. 





a us,A 




<r x io3 


OPE 


-0.144 ±0.024 


-0.028 ± 0.021 


-7.79 ±0.14 


Pheno. 


-0.135 ±0.003 


-0.028 ± 0.004 


-7.77 ±0.08 



spectral functions are very similar, but the errors are much lower using the phenomenological 

(k I) 

parametrization. Therefore, using phenomenology to describe the J=L component of SRr ' 
reduces considerably the theoretical uncertainty. 

We use as input value m s {2 G eV) = (9 5 ±20) MeV which includes the most recent determina- 
tions of m s from QCD Sum 

RulesEHEEEa 

lattice QCD ^ and r hadronic date 
With this strange quark mass input and taking and <5i?l°'°^ L from the QCD OPE and 

phenomenology respectively, as explained above, one can calculate 8R+ in @ from theory 

5R®& = (0.162 ± 0.013) + (6.1 ± 0.6) m 2 s - (7.8 ± 0.8) m 4 s = 0.218 ± 0.026, (11) 



H BO ID! III [EH 113 III 



where m s denotes the strange quark mass in MeV units, defined in the MS scheme at 2 GeV. 

In order to obtain a value of \V US \ from (|1U|) . we also need the experimentally measured 
Cabibbo-allowed R^y +ji and Cabibbo-suppressed R^'s contributions to the r decay rate. 

OPAL has recently updated the strange spectral function nJH In particular, they measure 
a larger branching fraction B{t~ — > K^tt + tt^i') which agrees with the previously one measured 
by CLEO©. From the OPAL data and using (JTU|) we get the result 

\V US \ = 0.2208 ± 0.0033 cxp ± 0.0009 th = 0.2208 ± 0.0034; (12) 

where we have used as input the PDG value for \V u d\ = 0.9738 ± 0.0005. The most important 
feature of this determination is the small theoretical error, that leaves the final uncertainty dom- 
inated by the experimental input. One can expect then that the better data sets from BABAR 
and BELLE will reduce significantly the error of our determination by improving the experi- 
mental input. Nevertheless, alread y now, our result is com petitive to the standard extraction of 
| V us | from semileptonic kaon de cay ^^^^^^^O l S ] anc j a new determination from fx / fn 
as extracted from lattice ^^ l 

One can use the value of |V^ S | thus obtained in (fT2*|) and determines the strange quark mass 
from higher moments. The weighted average of m s calculated for the different moments gives 

m s {M T ) = 84 ±23 MeV => m s (2 GeV) = 81 ± 22 MeV . (13) 

For details about the results for the different moments and the sources of the uncertainties se 

In our previous analysis based on the ALEPH data^J it was observed that m s displayed 
a strong dependence on the number of the moment k, decreasing with increasing k. With the 
recent CLEO and OPAL results finding a larger branching fraction B(r~ — > K~it +/ k~v), this 
dependence is much reduced, although still visible. This issue needs to be clarified with the help 
of better experimental data. 

4 Simultaneous fit of m s and |V^ S | 

The ultimate procedure to determine both \V US \ and m s from r hadronic decays, will be a 
simultaneous fit of both to a certain set of (k, I) moments. A detailed study including theoretical 
and experimental correlations will be presented elsewhere^l 



j n liO| we restrict ourselves to a simplified approach where all correlations were neglected. 
For this simultaneous fit of | V us j and m s we use the five OPAL moments^ from to R^' ^ . 

The central values we find from this exercise are 

\V US \= 0.2196 and m s (2 GeV) = 76 MeV . (14) 

These values are in very good agreement with our previous results in (|12j) and l)13j) . We expect 
that the uncertainties on these results will be smaller than the individual errors in (j!2j) and (j!3[) . 
but only slightly since the correlations between different moments are rather strong. 



5 Conclusions and remarks 

Using the strange spectral function updated by opalEI we get 

\V US \ = 0.2208 ± 0.0033 exp ± 0.0009 th , (15) 

and 

m s (2GeV) = 81 ± 22 MeV. (16) 

This result is expected to be highly improved in the near future due to the fact that the error 
is dominated by the experimental uncertainty and that uncertainty can be reduced with better 
data samples from BABAR and BELLE. But already now, the high precision r decay data 
from ALEPH and OPAL at LEP and CLEO at CESR provide competitive results for \V US \ 
and m s . The combined fit to determine both quantities including theoretical and experimental 
correlations is underway. 

The actual status of the \V US \ determinations has been nicely reviewed recently i n I33ISHl 
Though the CKM unitarity discrepancy has certainly decreased with the new theoretical and 
experimental advances, the situation is not yet as good as one could wish, and an accurate 
determination for \V US \ with the eventual precise measurement of the strange spectral function 
at BABAR and BELLE is desirable. With the value of |V^ S | in ((T2|) and using the PDG value 
| Kid I = 0.9738 ± 0.0005, one finds 

1 - \V ud \ 2 - \V US \ 2 - \V ub \ 2 = (2.9 ± 1.8) • 10~ 3 , (17) 

so that unitarity is violated only at the 1.6 a level. 

There are some open questions that will have also to be addressed'^. The (k, 0)-moment 
dependence of the m s prediction has been reduced after the recent OPAL and CLEO analyses 
finding larger branching fractions for t~ — > K~ir + ir~i>. Accurate experimental data will clarify 
the origin of the remaining moment dependence and will allow a consistency check of the whole 
analysis. 

Another point to be checked is the fulfilment of the quark-hadron duality between the 
QCD OPE and the OPAL spectral function. A first look at this question was presented in ^Sl. 
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